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Pion condensation in quark matter with finite baryon density 
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The phase structure of the Nambu - Jona-Lasinio model at zero temperature and in the presence 
of baryon- and isospin chemical potentials is investigated. It is shown that in the chiral limit and 
for a wide range of model parameters there exist two different phases with pion condensation. In 
the first, ordinary phase, the baryon density is zero and quarks are gapped particles. In the second, 
gapless pion condensation phase, there is no energy cost for creating only u-or both u and d quarks, 
and the density of baryons is nonzero. 
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According to a well-known point of view [lj, pionic degrees of freedom and, especially, the pion condensation 
phenomenon might play a significant role in the description of different nuclear matter effects, heavy-ion collisions, 
and the physics of compact stars. In reality, dense baryonic matter obeys an isospin asymmetry, i.e. the densities 
of up- and down quarks are different. In order that QCD adequately describes the isotopically asymmetric matter, 
usually the isospin chemical potential fii is introduced into the theory. It is well-known that in QCD with finite 
fij a nonzero pion condensate is generated if [ii is greater than the pion mass m n . This result was obtained in 
the framework of an effective chiral Lagrangians approach Q as well as in QCD lattice calculations, performed at 
zero or small values of the baryon chemical potential [is 01- However, these two approaches are not applicable for 
the description of an isotopically asymmetric matter at moderate baryon density. To overcome this problem, it was 
recently proposed to study the QCD phase diagram on the basis of Nambu - Jona-Lasinio (NJL)-type models, which 
contain quarks as microscopic degrees of freedom, in the presence of a baryon chemical potential \xb and an isospin 
fii one. In this way the influence of /is, /ij on both the chiral symmetry restoration effect and the formation of 
color superconducting (CSC) dense baryonic matter [f| was considered, but with a zero pion condensate (for earlier 
investigations of the NJL phase diagram in the plane of temperature-baryon density, but without isospin asymmetry, 
i.e. [ii = 0, see, e.g., [(|). In addition, it was found that if the influence of \ii is combined with the condition of 
electric charge neutrality, then essential changes of quasiparticle dispersion relations might occur in CSC matter. 
As a consequence, it is transformed into a so-called gapless CSC matter, since new gapless quasiparticle degrees of 
freedom appear (see, e.g., reviews Q). Recently, the pionic condensation effect was investigated in some NJL models 
at nonzero values of fiB and 8, 9]. However, in these papers the existence of a homogeneous dense baryonic phase 
with pion condensation was not investigated. So, the question whether it is possible to describe such a phase in the 
framework of NJL models is opened up to now. 

To answer this question, we study in present paper the phase structure of the NJL model in terms of \ib and 
We show that, depending on the parameters of the chosen NJL model, a homogeneous phase of an isospin asymmetric 
quark matter with pion condensation and with finite baryon density may appear. Since quarks are gapless in this 
phase, we call it gapless pion condensed (GPC) phase. In addition, a pion condensed phase with zero baryon charge 
is also presented in the phase structure of the model. In contrast to GPC phase, quarks are gapped in it. In the 
following we restrict ourselves to values /is < 1 GeV in order to avoid the inclusion of a Cooper pairing of quarks 
into consideration. (At greater hb the diquark condensation should be taken into account |lfJ.lll|.') 

Let us study the properties of quark matter at moderate baryon and isospin densities and at zero temperature in 
the framework of the NJL model with two quark flavors 

L q = q 1 v id v - m + /!7° + 5flT 3 j° q + G {qqf + {qrfrqf , (1) 

where /i = /is/3 is the quark number chemical potential and the quantity S/i is a half of the isospin chemical potential, 
i.e. S/j. = fii/2. In l[T|l n (i = 1, 2, 3) are Pauli matrices. Evidently, at nonvanishing current quark masses, m ^ 0, the 
Lagrangian Q is invariant under the baryon Ub(1) symmetry and the parity transformation P. Moreover, without 
the /ir-term it is also invariant under the isospin SUi(2) group which is reduced to Ui 3 (l) at fj,j ^ 0. Usually, the 
properties of the ground state of the model are characterized by two order parameters. The first one is the quark 
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FIG. 1: The phase portrait of the model at m = in terms of /x, Sfj, for the set 1 of parameters, A = 0.65 GeV, G = 5.01 GeV - . 
The region D corresponds to the fully symmetric phase with M — 0, A = 0. A, B, and C are phases with pion condensation. 
The solid lines fii c (Sfi) and ^2c(SfJ,) are the critical curves of first order phase transitions. The dotted line, A (5/j,), is the pion 
condensate at [i = 0. The baryon density is zero in the phase B and nonzero in A, C, D. There are gapless quarks in the phases 
A, C, D. In the phase B quarks are gapped. 



condensate (qq) which is responsible for the chiral symmetry breaking and does not spoil the parity and isotopical 
symmetry. The second one, called pion condensate, is the ground state expectation value of the form (<?7 Tig). If 
it is nonzero, spontaneous breaking of parity and isotopic symmetry occur in the system. 1 At /!/ = only the 
quark condensate is relevant to the physics of the model, but at nonzero /i/ it is desirable to take into account both 
condensates. In the last case the competition between them is governed by the thermodynamic potential (TDP) which 
in the mean field approximation has the following form which can be obtained, using, e.g., the technique of 

(M~m) 2 + A 2 v-^ f d 3 p , 
A) = < JL ^J^\E a \ = 

(M - to) 2 + A 2 f d 3 p 



P)wl 

m 



The summation in the first line of J2J runs over all quasiparticlcs (a = u, d, u, d), where 

E u = E A — /i, E a = E A + fi, 

Ed = E A - M , Eg = E A + /x, (3) 

and E A = y/ (£' ± ) 2 + A 2 , E ± = E ±Sfi and E = yjp 2 + M' 2 . The factor 3 in front of the summation symbol in 
(0) indicates the three-fold degeneracy of each quasiparticle in color. Moreover, in order to avoid usual ultraviolet 
divergences, the integration region in @ is restricted by a cutoff A, i.e. \p\ < A, where A < 1 GeV. The gap 
coordinates of the global minimum point of the function f2(M, A) are connected with condensates in the following 
way: 

M = to - 2G(qq), A = -2G l (g 7 5 r l9 ), (4) 



1 One may, in addition, suppose that (<77 T 2<?) is also nonzero. Then, using the Ui 3 (l)-rotations, this quantity can be reduced to zero. 
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FIG. 2: The gapless phase A of the NJL model at m = 0. The notations are the same as in Fig. 1. The schematic curve L 
divides A into the region with only u gapless quarks (to the right of L) and the region with both u and d gapless quarks (to 
the left of L). In the triple point a three phases A, B, D coexist. 



so if A is nonzero in the global minimum point (GMP), then parity as well as isotopic symmetry are spontaneously 
broken and the pion condensation phase is realized. Note that the quark gap M is just the dynamical quark mass. 
From @ it is possible to obtain the gap equations 

Q _ 0n(M,A) ^ M-m 6M f d 3 p ( 8(E+-n)E+ , 0(E A — fx)E~ 



dM ~ 2G J (2tt) 3 £1 E+ E a 



dCl(M,A) _ A 6A r d 3 p ( 0(E+-») | g(£ A -/i) j (5) 



dA 2G J (2tt) 3 I E+ E A 



We will investigate the phase structure of the NJL model both at /i/ ^ and \x ^ 0. If the current quark mass m 
is not zero, then the pion condensation in the framework of NJL approach is occured at /i/ > 0,0. Just for these 
values of fii it is reasonable to search for the pion condensed phase with finite baryon density. In general, it is a rather 
hard task, so in the present paper we restrict ourselves by the zero current quark mass case only (m = 0, irij = 0), 
thus studying the phase structure of the model for all fii > = m T . At m — and fii ^ the symmetry group of 
the initial model is Ub(1) x Ut 3 l(1) x Ut 3 r.(1) x P. The system of the gap equations (J5J has three different solutions: 
i) M = 0, A = 0, ii) M ^ 0, A = 0, iii) M = 0, A / 0. Having all the solutions of the gap equations, it is necessary 
to compare the values of the TDP on them, selecting by this way the minimal one, i.e. the GMP. If the GMP of 
the TDP is achieved on the solution i) , then the ground state of our system is invariant under the initial symmetry 
group of the model. The solution of the type ii) corresponds to the phase, symmetrical under Ub(1) X Ui 3 (l) x P 
group, where Ui 3 (l) acts in the following way: q — > exp(iaT3)q. Finally, if the GMP of the TDP lies in the point of 
the form iii), then in the ground state the quark condensate is zero, but the pion condensate is nonzero. In this case 
the pion condensation phase is realized in the model, and the ground state is symmetrical under Ub(1) X Uai 3 (1), 
which transforms left- and right quark flavor doublets in the following way: qh — > exp(iar3)gi, qji — > exp(— iar^qn. 
Generally, at m — 0, due to a symmetry of the TDP, it is sufficient to study its GMP only in the region M, A > and 
for nonnegative values of Sfi, fi. For the moment let us fix the parameters as G = 5.01 GeV -2 , A = 0.65 GeV (set 1). 
The dependence of quark matter properties on the different sets of model parameters G, A will be discussed later. 

Following the way, described above, we obtain for set 1 of parameters the model phase structure that is represented 
in Fig. 1 in terms of Sfi = fii/2 > and \i. There, the region D corresponds to the totally symmetric phase of the 
model, in which M = and A = 0. This phase is separated from the regions A, B, C by the critical curve /i = /ii c ((5/i) 
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FIG. 3: Typical behaviour of quark energies in the phase A and to the right of L (see Fig. 2). The case presented corresponds 
to <5/i = 0.06 GeV, fi ~ 0.306 GeV. In this case the pion condensate takes the value A(/i, S/i) ~ 0.303 GeV. 

of first order phase transitions. In A, B, and C the GMP of the TDP has the form iii), so these regions correspond 
to phases with pion condensation. 

In the phase B the pion condensate depends only on <5/x (not on fj,), so A = A (S/i), where A (i5/i) is the pion 
condensate at /i = 0. Since for all points of the region B we have the relation /.i < A (<5/i) (see Fig. 1, where A (8fi) 
is depicted by the dotted line), it is clear that the value of the TDP in the corresponding GMP point, f2 m i n , does not 
depend on the quark chemical potential n (this follows from the second line of J2J). So, in the phase B the baryon 
density ns is identically zero (jib = — <9f2 m i n /(3<9/i)). In contrast, in the phase A, which is presented in more details 
in Fig. 2, the pion condensate depends both on fi and S[i, i.e. A = A(/x, Sfx). Hence in the phase A the baryon density 
is nonzero. Since on the curve /i = /Z2 C (^A*) of Figs 1,2 we have A(/i, Sfi) ^ A G (5/x), the transition from the region 
A to B is a discontinuous first order phase transition. In the so-called triple point a of Fig. 2 three phases A, B, D 
coexist. In the phase C that is separated from B by the dotted line [i — A (5fi), the pion condensate also does not 
depend on fi. Moreover, here it is the same as in the phase B, i.e. A = A (5)i). However, since for the points of the 
region C we have fi > A (5fi) (see Fig. 1), the minimal value S^i min of the TDP @ is changed with \i. As a result, 
the baryon density in the phase C is nonzero as well. Note that the phase C is located at values <5/i > A. In this case 
one might believe that it is an artefact of the NJL model and not adequate to the real situation. On the other hand, 
the phase A corresponds to rather small values of /i and /i/, so it might be relevant to real physics and needs some 
more discussions. 

As it is easily seen from Fig. 2, for the set 1 of G and A the height of the region A, i.e. the quantity H = /i lc (0)—A o (0) 
(in the model under consideration, at m = the value A o (0) is nothing else than the dynamical (constituent) quark 
mass M q m 0.3 GeV at /i = 0, £/i = 0), is approximately 15 MeV, whereas the extention of it along the 5/x-axis 
reaches 70 MeV (for /i/ = 140 MeV). Apart from the baryon density, there is another physical characteristic which is 
different in A and B phases. It is just the minimal energy, necessary for creating a single quark in the ground state. 
Since the quark condensate is equal to zero both in A and B phases, one should take M = in formulae 10 in order 
to obtain the energies of u- and d quarks in the pion condensed phases A and B. So, 

E u = VCp - + A 2 - n, E d = v/(p + ^) 2 + A 2 - /i, (6) 

where p = \p\. If A > /i, the minimal values of the quantities E u ^ are greater than zero. This fact means that for 
creating a single quark in the phase B, where A = A Q ((5^) > /i, one needs an energy that is greater than some finite 
(nonzero) amount (also called gap). Hence, quarks are gapped in the B phase. In contrast, for all points of the phase 
A the relation A = A(/i, Sfi) < [i is fulfilled. Hence, in this phase the minimal energy necessary for creating u- or 
both it- and d quarks is always zero. So, in the phase A we have gapless quarks. The line L that is defined implicitly 
by the equation /i = \J (<5/i) 2 + A 2 (/i, 5ji) divides the region A into two parts (in Fig. 2 the line L is represented very 
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FIG. 4: Energies of u quarks vs three momentum p for three different points of the phase A that are to the left of the curve L 
(see Fig. 2). I - the case of Sfi = 0.01 GeV, y w 0.31 GeV, A(n,Sfj,) « 0.298 GeV. II - the case of 5fx = 0.03 GeV, fi « 0.31 
GeV, A(/i, <5^) « 0.297 GeV. Ill - t 
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FIG. 5: The analogue curves as in the Fig. 4, but for the d quarks. 

schematically). To the right of L only u quarks are gapless in the phase A. Typical behaviours of quark energies vs 
p in this region are shown in Fig. 3. (Note, \E U \ turns into zero at two different values of p.) To the left of L both 
u- and d quarks are gapless ones. This situation is illustrated in Figs 4 and 5, where for /i = 0.31 GeV the energy 
of the u quark is depicted for three different values of Sfi (Fig. 4), whereas in Fig. 5 the energy of the d quark is 
represented for the same values of /i and <5/i. It is easily seen that to the left of L the energies of u- and d quarks turn 
into zero only in a single point of the p variable. Due to these features, it is naturally to call the phase A the gapless 
pion condensed (GPC) phase. (Note, in the phases C and D there are also gapless quasiparticle excitations.) 

Up to now in considering the GPC phase we have dealed only with the set 1 of parameters G and A. For this set 
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FIG. 6: The heights H of the gapless phase A vs cutoff A for different fixed values of constituent quark masses M q (in GeV). 

the constituent quark mass M q in the vacuum, i.e. at zero /i and /ij, is approximately 0.3 GeV. It follows from the 
relation, connecting G, A, and M q |T3| : 



Using (J2J, it is possible to study the dependence of the height H = /ii c (0) — A o (0) of the GPC phase A on the 
parameters A and M q (see Fig. 6). Clealy, positive values of H indicates on the presence of GPC phase. In this 
figure one can see a set of curves. Each curve represents the behaviour of H vs A at fixed M q . (The result H = 15 
MeV, obtained above for the set 1, corresponds to the curve with M q = 0.3 GeV at A = 0.65 GeV.) It follows from 
Fig. 6 that the GPC phase is presented in the phase structure of the NJL model for a rather wide range of model 
parameters. The smaller M q , the greater the height of the GPC phase at fixed A. In this context, it is worth noting 
that some parametrizations of the NJL model correspond to rather small values of M q . For example, in [14] the 
constituent quark mass was taken as small as 0.233 GeV, whereas A = 1 GeV. Hence, the height of the GPC phase 
in this parametrization reaches the value 50 MeV. In contrast, in the paper the quantities M q = 0.428 GeV and 
A = 0.58 GeV were used for the consideration of the NJL model at nonzero /ib amd fii. It is clear from Fig. 6 
that for such values of parameters the GPC phase is absent in the phase structure, so it is no wonder why in ||| the 
homogeneous phase, describing quark matter with finite baryon density and pion condensation, was not found. 

In summary, we have investigated the phase structure of the two flavor NJL model at finite baryon and isospin 
chemical potentials for a current quark mass m taken to be zero. It was shown that at the cutoff value A = 0.65 
GeV and for the constituent quark mass M q = 0.3 GeV (or G = 5.01 GeV -2 ) there are three different phases in the 
model, A, B, D, at < 5fx < A. The phase D corresponds to quark matter with finite baryon density and zero pion 
condensation. In the phases A and B pions are condensed. However, baryon density in the phase B is zero, but in 
A it is nonzero. Moreover, the phase B is a gapped one, since creation of a single quark in the ground state of this 
phase necessarily costs some finite amount of energy. In contrast, in the phase A u-quarks or both u- and (i-quarks 
are gapless. These three phases can coexist in the triple point a (see Fig. 2), since transitions from one phase to 
another are discontinuous. 

It is well-known that pion condensation takes place if /xj is greater than the pion mass m n (see, e.g., |9j). For 
simplicity, in the present paper we have considered the zero current quark mass case, where pion mass is zero. That 
is why the GPC phase is observed at rather small values of fij. In a more realistic situation with m ^ the pion mass 
equals approximately 140 MeV. We believe that in this case GPC phase might be predicted as well, but at fij > m n . 

In conclusion, there arises the interesting question, whether a pion condensate with finite baryon number density 
might exist in nature. 
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